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Abstract 

"Oscillations" occur in quite different kinds of many-particle-systems when 
two groups of particles with different directions of motion meet or intersect at a 
certain spot. In this work a model of pedestrian motion is presented that is able to 
reproduce oscillations with different characteristics. The Wald-Wolfowitz test and 
Gillis' correlated random walk are shown to include observables that can be used 
to characterize different kinds of oscillations. 

1 Introduction 

The widely used word oscillation is defined for this article in a very general sense 
as the phenomenon that the direction of flow changes in certain intervals when two 
groups of particles compete for the chance to move at a narrow bottleneck. Oscillations 
in pedestrian dynamics have been dealt with before (see [ 1 1 for example). The aim of 
this work however is to give a framework for an analytical and quantitative treatment of 
oscillations in observation, experiment and simulation, which has not been done before. 
There are three extreme types of oscillation: 1) the zipper principle: the flow direction 
changes after each particle. 2) no oscillations: the flow direction only changes after 
one group has completely passed through the bottleneck. 3) uncorrected oscillation: 
the statistics of the flow direction is the same as for coin-tossing experiments. While in 
the first two scenarios the whole process is completely determined as soon as the first 
particle has passed the bottleneck, there is absolutely no influence from one event to 
the next in the third scenario. Between those extremes there is a continuous spectrum 
of correlated oscillations with an influence from one event to the other but no full 
determination. The outline of the paper is as follows: First we present the parts of our 
model of pedestrian motion that are crucial for oscillations. We then describe the test 
scenario in detail. The third step is the presentation of the observables that characterize 
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a discrete oscillation. Finally we present the simulation results for different sets of 
parameters. 

2 The model 

A model of pedestrian motion is presented now that expands the model of [3 1 for 
velocities v > 1 and adds some more features. Here we will only present the elements 
of the model that are involved in the phenomenon of oscillations. The full model in all 
details will be presented elsewhere [4 |. The model is discrete in space (cells) as well 
as in time (rounds). Space is discretized into a regular lattice with cells with border 
length 40 cm. Each cell may be occupied by at most one agent as we will call the 
representation of a real person in the model. At the beginning of the simulation the 
distance of each cell to each exit is calculated and stored in the integer valued static 
floor field S^ y . So S% y is the distance of cell (x, y) to exit e. Additionally to the static 
floor field there is an also integer valued dynamic floor field D xy = (D xy , D^y) which 
in our model is a vector field. Whenever an agent leaves cell (x, y) and moves to cell 
(x+i, y+j) the vector (i, j) is added to the dynamic floor field at the position (x, y). In 
between the rounds of the agents' motion the dynamic floor field is subject to diffusion 
and decay with the two components diffusing and decaying independently: each entry 
(positive or negative) in a component decays with probability 5 and if it does not decay, 
it diffuses with probability a to one of the neighboring cells. 

The agents can have different maximal velocities. This is implemented by the number 
of cells an agent is allowed to move into any direction during one round. As one round 
is interpreted as one second, typical velocities are two to five cells per round. See 1 5 1 
for details. So for each agent a there is a set R a of cells, which are unoccupied and 
visible (not hidden by a wall or another agent) and which can be reached during the next 
round. At the beginning of each round each agent a randomly chooses a destination 
cell (x°i, Vd) out °f R a - The probability for a cell (x, y) within R a to be chosen as 
destination cell by an agent currently positioned at (xq, y%) is: 

pl <y ex exp (k a s S e = f[a ^ y ^ x '' y ' ] + k a D ((x - x a )D^ y . + (y- %°)^g >2/ g)) (D 

Compared to the model variants in |6| where desired trajectories are calculated in ad- 
vance, by applying the same update rules v max times and then different possible ways 
to use that trajectory to reach the final cell are investigated, here only a desired destina- 
tion cell is calculated and the actual movement process is not implemented as planning 
process but as "real" motion. Note especially the equality of all cells of R a , where in 
principal the relative position of a cell to an agent does not matter. 

The dynamic floor field therefore influences the motion via the scalar product with 
a possible direction of motion. g^^^-^o'Vo^^^^) means mat m me p resence f more 
than one exit the choice of an exit is a function of the number (group affiliation) and 
the position of the agent as well as the position of the exit. Here all agents are only 
allowed to use the exit on the other side of the bottleneck compared to their starting 
position. 



After all agents have chosen their destination cell in parallel, the agents start to 
move sequentially. It is important to mention that on their way to their destination cell 
they execute single steps within the Moore neighborhood [7| which implies that for 
v = 1 one step within the Moore neighborhood is possible. The sequence of all single 
steps of all agents is fixed randomly at the beginning of the motion part of a round. 
Compared to the model treated in 1 8 1 there is an additional shuffle as an agent does not 
necessarily execute all his steps at once. Whenever an agent has stepped on a cell and 
moved on, the cell remains blocked for all other agents until the rest of the round. This 
represents the speed-dependent occupation of space during motion |9lQ0||6), which is 
essential for a realistic fundamental diagram. Within a Moore neighborhood an agent 
always chooses the cell which lies closest to his destination cell and which has not been 
blocked before by another agent. If there is no unblocked cell left which is closer to the 
destination cell than the current cell, the round ends for that particular agent. Therefore 
an agent can fail to reach his destination cell. If two or more agents choose the same 
cell as destination cell or need the same cell as part of their path to their destination 
cell in the calculations for this work simply the order of movement decides about the 
winner of such a "conflict". 



3 The scenario 
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Figure 1: Initial positions of the agents in the scenario (schematic view). 



At the beginning of the simulation two groups (A and B) with Na = Nb = 202 
(N = Na + Nb) agents each are placed in two rows as shown in figure Q. They have 
to pass the bottleneck with a width of one cell in the middle and proceed to the exit on 
the other side. An agent is counted as having passed as soon as he leaves the bottleneck 
cell in the direction of his exit. The possibility that he moves back into the bottleneck is 
given but neglected. Only the first n = 100 agents passing the bottleneck are taken into 
account to make sure that there are always two almost equally sized groups competing 
for passage. Of these n = 100 agents the number of agents of group A is called riA 
and of group B ng. n can not be chosen arbitrarily large, as at some size the crowds 
in front of the bottleneck become that dense that no passing at all is possible. 

There are no new agents entering the scenario. 



4 Observables 



4.1 The Wald-Wolfowitz test 

The Wald-Wolfowitz test can be used to check for correlations in time series whenever 
the statistic is dichotomous or can be dichotomized, i.e. when two kinds of events or 
data are present. The probabilities for the events to happen do not need to be equal. In 
our case those two kinds of events are "A member of group A passes the bottleneck" 
and "A member of group B passes the bottleneck". The Wald-Wolfowitz test makes 
a statement about the expectation value and the variance of the number of runs (and 
other observables) if there are no correlations between the different events. In the 
terminology of the test a run is a series of events of one kind not interrupted by an 
event of the other kind. 

The crucial quantities of the Wald-Wolfowitz test are: 

• (R(n Al tib)) is - for given ua and tib - the expectation value of the number of 
runs for an uncorrelated series of - in our case - bottleneck passages. 

• an is the standard deviation of the number of runs. 

• r is the number of runs in a certain simulation, 

• z = r ~ W"-4,ra_B)) j s ^ standardized test variable. 

or 

For nA,ns S> 1 (at least nA,ns > 10) the distribution of the number of runs 
becomes comparable to the normal distribution which implies that the null hypoth- 
esis "The events are uncorrelated" can be rejected on an a-level of significance if 
\z\ > q(l — a/2), with q(l — a/2) as quantile of the standard normal distribution for 
the probability 1 — a/2. For smaller nA,ns one has to make use of tables to decide 
about the rejection of the null hypothesis 1 1 11 1121 . 



For an uncorrelated process one has for given ua, tib 

It>( W n n An B , 1 n A (n-n A ) 

{R{n A ,n B )) = 2 h 1 = 2 hi (2) 

n n 



2n A n B {2n A n B - n) = ((fl) - - 2) 

n 2 (n — 1) n — 1 



(3) 



The Wald-Wolfowitz test is independent of underlying distributions; therefore it 
does not make use of deviations in the distribution of n A in subsequently repeated 
simulations of a scenario. This additional information can be made use of by mapping 
the process onto a correlated random walk. 



4.2 Correlated random walk 

In correlated random walk models 1131 1141 the probabilities for the direction of the 
next step depend on the direction of the last step. The "correlated random walker" 
keeps his direction of motion with probability p and changes it with probability 1 — p. 
In the first step the direction is chosen with equal probability. In the thermodynamic 
limit (n — * oo) one then gets a normal distribution Af(n/2, jz—& A (p = 0.5)) for the 



probability that the walker made ua of n steps to the right 1 1 5 j . Here M(Xq , a 2 ) is the 
normal distribution with maximum at Xq and variance a 2 , and a 2 A (p = 0.5) = n/4 is 
the variance for the number n A of steps to the right in the case of uncorrelated random 
walk (p = 0.5). The position x n of the walker after n steps is ua — tib- Numerical 
calculations (see figure|3]in appendix|X|i show that for n = 100 the thermodynamic 
limit is a good approximation (relative error for the standard deviation < 1%) for 0.2 < 
p < 0.8 The probability p to keep the direction should equal the sum of the correlation 
coefficients to find an event A directly followed by an event A and B directly followed 
by B. 

P ~ CCkeep direction = CCfe = P({AA) Aevents=l) + P {{BB) Aevents=l ) (4) 

We map the events of the oscillation experiment on a special correlated walk - 
Gillis' random walk in one dimension 1 16l : 

• "An agent of group A passes the bottleneck." — > "Random walker moves one 
step to the right." 

• "An agent of group B passes the bottleneck." — » "Random walker moves one 
step to the left." 

Because of equation (|2ji there is a connection between (R(ua, ns)) - the expecta- 
tion value of the number of runs - and a a- Since a a — <?b one has 



5 Values for an uncorrelated process 

For n — 100 and Na = Nb = 202 the expectation value for ua has to be 50 (= n/2) 
due to the symmetry of the scenario. Even if a finite group size effect would in a strict 
sense imply a correlation, we here want to give not only the values for an uncorrelated 
process with no finite group size effect, but also for an else uncorrelated process with 
finite group size effect. In the first case, where the passing order can be imagined to 
be determined by a guardian at the bottleneck who decides about the right of passage 
of some member of one of the two groups by throwing a coin, the basic distribution is 
binomial. (All numerical values are calculated for n — 100 and Na = Nb = 202) 



(*( ffA ))«-(n +ffA )(2_ ffA ) 
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(R) = P(n A ){R(n A , n B )) + 50.5 



(8) 



The latter case, where the guardian decides about the right of passage with equal prob- 
ability for each of the remaining individuals, is governed by the hypergeometric distri- 



bution. 
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4.34 



50.62 



(9) 

(10) 
(11) 



It is assumed that the size of groups on both sides of the bottleneck has no influence 
on the movement probability as long as both groups have some minimal size which 
should be exceeded by construction during the counting process as the largest possible 
inequality is 102:202. However even if this was not the case, this comparison shows 
that the expectation values of the two possible cases of an uncorrected process differ 
only slightly. We now can take these results as basis for comparison to the results of 
our simulations. 



6 Oscillations without dynamic floor field 
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Figure 2: One of twelve possible situations after an agent has passed. 



For kjj = there is some sort of "intrinsic" correlation, which stems exclusively 
from the hard-core exclusion of the agents. As the last agent who passed, occupies one 
of the three cells which are adjacent to the bottleneck cell, there are on average fewer 
agents on the adjacent cells on the side of the group that could not move than on the 
other side. Figure[2]shows the situation when there is the maximum of three agents of 
the group that moved last and the maximum of only two agents of the other group. Of 
course there can be other combinations (3:2, 3:1, 3:0, 2:2, 2:1, 2:0, 1:2, 1:1, 1:0, 0:2, 
0: 1 , 0:0) and the group that moved last can have fewer agents adjacent to the bottleneck, 
but on average there will be more. Therefore on that side there will on average also 
be more agents that plan to move to the bottleneck cell and so also the probability that 
one of them wins a possible conflict is increased, as the winner of a conflict is chosen 
with equal probability out of all competitiors. Note that this effect crucially depends 



on the choice of von Neumann neighborhood or Moore neighborhood as neighborhood 

forv max = 1. 

For ks = 0.5 we repeated the calculations for the three cases v max = 1, v max = 3 
and v max = 5 one-hundred times each. The tables of this and the following sec- 
tion show the average number of runs that occured, the average number of runs the 
Wald-Wolfowitz test predicts for uncorrelated behavior, average z-values as well as 
the number of times (out of 100) the z-value was smaller than -1.95 (5%-significance 
level) and smaller than -2.58 (1%-significance level), the standard deviation of tia, the 
evacuation time (e.t.) in rounds and the probability to continue in the same direction as 
calculated from a a (p) as well as calculated from correlation coefficients (cck). 
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Note that (R(ua, ns)) is the average of the values for (R(tia, calculated from 
the 100 simulation results for ua respectively using equation 10 and averaged in 
the way it is done in equation (|8j however using the distribution of simulation results 
and not a theoretical distribution P(ua)- The other averages are also averages over 
100 simulations. 



7 Oscillations with dynamic floor field 

For all simulations ko = 0.3 and ks = 0.5 has been set. This is quite a small value, 
however for ks = 1 .0 or even larger the crowds on both sides of the bottleneck become 
too dense, such that it becomes difficult for an agent to pass through them after he has 
passed the bottleneck. For each set of parameters the simulation was repeated 100 
times. We mainly experimented with the maximum speed v max (which was equal for 
all agents in a scenario) and the parameters for diffusion a and decay 5. The largest 
deviations from the uncorrelated case are marked in the following tables. 
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(Also see figures|4][5][6][Zl|8]and[9]in appendixlAh Other simulations showed that for 
8 > 0.1 the effect begins to vanish as the dynamic floor field decays too fast to have 
a significant influence and for 5 < 0.01 the effect is hidden by the dominant dynamic 
floor field that then induces irrational behavior. The agents begin to move in circles 
instead of heading to the bottleneck. Note that effects like this can also occur for too 
small a, if the dynamic floor field becomes too strong on some cells. 

Compared to the case without dynamic floor field, the dynamic floor field with 
its time dependence brings in an another dimension that can be analyzed: As long 
as the dynamic floor field has not reached a steady-state, the oscillation will be time 
dependent. This implies that the measured values either change with n or if n is kept 
constant but the measurement process starts not with the first event but later. For a — 
1.0, v max = 1 and n = 150 we got z = —3.92, which is quite a difference to z = 
—3.26 for n = 100. However except for the case a = 0, where they became slightly 
larger, the correlation coefficients did not change at all. As the existence of time- 
dependent oscillation types is model-dependent, the analysis of this phenomenon shall 
not be done in further detail in this work. 

8 Discussion 

Some observations made in the ku = 0.3 data are as follows. 

• The evacuation time (in rounds) drops dramatically compared to kn = simu- 
lations. So it can be assumed that the larger groups of agents passing the door 
represent a more efficient behavior. This can be interpreted as less time being 
consumed by conflict solution processes at the bottleneck. 

• The minimal z-value for each v max is found at larger a for larger v max . (See 
also figure|4]) 

• The same holds for the minimum of the evacuation time. (See also figure|5]) 

• Correlations manifest themselves for v max = 1 typically in smaller z-values and 
for v max — 5 typically in larger <ja- 

• a a for v max = 3 and v max = 5 are that large that (nA,riB > 10) is not always 
fulfilled. 

• For simulations with dynamic floor field, the correlations are stronger for v max = 
3 and v max — 5 than for v max — 1. This is due to the larger area of influence 
for higher speeds: more agents could choose the bottleneck cell as their desti- 
nation cell. But due to the direction of the dynamic floor field mainly agents 
of the group that moved last indeed do choose the bottleneck cell as destination 



cell, which leads to a stronger outnumbering of agents attempting to follow one 
of their group on the bottleneck cell compared to agents trying to change the 
direction of the flow, than in the case of v max = 1. 

• Also for simulations with dynamic floor field p{cja) is most of the time much 
larger than cc k , however within each v max -set of results there appears to be a 
tendency that cc k and p are positively correlated. The difference points to a 
dynamic floor field that is too strong to reverse direction with only one agent 
passing. If for example after a sequence of agents of group A one agent of group 
B passes the bottleneck, a typical sequence will look like AAAAAABAAAA 
if the dynamic floor field does not change direction with that single agent. If 
however it does change, a typical sequence would be AAAAAABBBBB. In 
the first case a a will grow, as the dominance of agents of one group outlasts 
the accidental event with small probability that one agent of the other group 
passes. In the latter case a a will be comparatively small since one long run of 
As can be followed by an equally long run of Bs and vice versa. The correlation 
coefficient cc k however is much less affected by this phenomenon as it is not 
distinguished between AA and BB sequences. Take for illustration an A-B- 
symmetric n = 11 example: {AAAAAABAAAA, BBBBBBABBBB) has 
a A = 4.5 and cc k = 0.8 and (AAAAAABBBBB, BBBBBBAAAAA) has 
a a = 0.5 and cc k = 0.9. So while a a becomes larger with a stronger dynamic 
floor field, cc k becomes larger. Consequently did a simulation with v = 3, 
a = 1 and S = 0.1 result in cc k = 0.64 and p = 0.61 instead of cc k = 0.72 
and p = 0.95 as for the S = 0.03 simulation. This shows that a comparison of 
p and ccfe can give a hint that S was chosen too small, as for real pedestrians the 
passing of one individual is enough to completely change the odds. 

• Simulations without dynamic floor field gave values for a a and p that might 
make p look dependent of v. However the reason for this is a relatively broad 
distribution of the a a- In 20 additional v max = 1 calculations of a a with n = 
100 simulation repetitions each a a variied between a a — 5.34 and a a = 6.72 
this implies values for p between p — 0.53 and p = 0.64 with an average of 
p = 0.575. At the same time cc k remained relatively constant between cc k = 
0.62 and cc k = 0.64. That p even at average remains smaller than cc k probably 
results from a reduced local density directly in front of the bottleneck if a group 
had a sequence of agents passing the bottleneck. The reduced local density then 
increases the probability of a change of the flow direction. This is an effect of 
the kind of a hypergeometric distribution, that reduces the standard deviation 
(see above). Here too p is much more affected than cc k as cc k is not affected by 
the relative total numbers of ^4s and Bs within the first n events, but only by the 
number of changes of the flow direction. 

The significant drop in the evacuation time compared to k D = is a sign that the 
dynamic floor field not only leads to a more efficient motion into the bottleneck but 
also out of the bottleneck and through the group on the other side. 

It is not very surprising, that the effects of the dynamic floor field are strongest for 
larger v max at larger a, since the larger neighborhood in which faster agents can move 



during one round makes it necessary that also the dynamic floor field diffuses faster to 
keep up with the agents. 

If a a becomes too large, moderate z-values are no indicator that the events are 
uncorrelated. Contrary to that the scenario is just too small for such sets of parameters 
to make large z-values possible. 

9 Summary, conclusions and outlook 

A model of pedestrian motion was presented, which is able to modulate the character- 
istic observables of oscillations. The model was tested concerning this capability using 
the Wald-Wolfowitz test, a comparison to correlated random walks and a calculation of 
correlation coefficients. This framework for characterizing should be applicable to any 
model of pedestrian motion as well as corresponding experimental and empirical data. 
However these observables would show details that are far from having been measured 
in reality at present. Anyway the scenario of oscillations at bottlenecks did show up a 
rich variety of phenomena and observables, which might be used to characterize and 
distinguish models of pedestrian motion. The investigation of time-dependent oscilla- 
tion types as well as the distribution of the length of runs could be further elements in 
such discussions. Observations and experiments in the field would be of great interest. 
Since our model is able to reproduce quite different types of oscillations it probably 
will also be able to reproduce the range of oscillations that occur in reality. 
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A Figures 




Figure 3: Relative difference of the numerically calculated standard deviation for the 
number of steps into a certain direction for n = 100 steps against the theoretical ther- 
modynamical (n — > oo) limit of the same standard deviation in dependence of the 
probability p to continue motion in the same direction. 




Figure 5: Average evacuation times for different a at S — 0.03. 
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Figure 7: Correlation coefficients for different a and event distance at v max = 1 




Figure 8: Correlation coefficients for different a and event distance at v max — 3 
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Figure 9: Correlation coefficients for different a and event distance at v„ 



